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The measurements of the angular power spectrum of the Cosmic Microwave Background (CMB) 
anisotropy has proved crucial to the emergence of cosmology as a precision science in recent years. 
In this remarkable data rich period, the limitations to precision now arise from the the inability to 
account for finer systematic effects in data analysis. The non-circularity of the experimental beam 
has become progressively important as CMB experiments strive to attain higher angular resolution 
and sensitivity. We present an analytic framework for studying the leading order effects of a non- 
circular beam on the CMB power spectrum estimation. We consider a non-circular beam of fixed 
shape but variable orientation. We compute the bias in the pseudo-C; power spectrum estimator 
and then construct an unbiased estimator using the bias matrix. The covariance matrix of the 
unbiased estimator is computed for smooth, non-circular beams. Quantitative results are shown 
for CMB maps made by a hypothetical experiment with a non-circular beam comparable to our 
fits to the WMAP beam maps described in the appendix and uses a toy scan strategy. We find 
that significant effects on CMB power spectrum can arise due to non-circular beam on multipoles 
comparable to, and beyond, the inverse average beam-width where the pseudo-C; approach may be 
the method of choice due to computational limitations of analyzing the large datasets from current 
and near future CMB experiments. 

PACS numbers: 98.70.Vc,95.75.Pq,98.80Es 

I. INTRODUCTION 

A golden decade of measurements of the cosmic microwave background anisotropy has ushered in an era 
of precision cosmology. The theory of primary CMB anisotropy is well developed and the past decade has 
seen a veritable flood of data [1, 2]. Increasingly sensitive, high resolution, 'full' sky measurements from 
space missions, such as, the ongoing Wilkinson Microwave Anisotropy Probe (WMAP) and, the upcoming 
Planck surveyor pose a stiff challenge for current analysis techniques to realize the full potential of precise 
determination of cosmological parameters. As experiments improve in sensitivity, the inadequacy in modeling 
the observational reality start to limit the returns from these experiments. 

A Gaussian model of CMB anisotropy AT(q) is completely specified by its angular two-point correlation 
function. In standard cosmology, CMB anisotropy is expected to be statistically isotropic. In spherical 
harmonic space, where AT(q) = J2i m a lmXi m {c[), this translates to a diagonal (ai m a*, m ,) — Ci6u>5 mm > where 
Ci, the widely used angular power spectrum of CMB anisotropy, is a complete description of a Gaussian 
CMB anisotropy. Observationally, the angular power spectrum being a simple, robust point statistics is the 
obvious first target for cosmological observations. Theoretically, the C; are deemed all important since the 
simplest inflation models predict a Gaussian CMB anisotropy. In this case, the power spectrum provides an 
economical description of the CMB anisotropy allowing easy comparison to observations. 

Accurate estimation of Ci is arguably the foremost concern of most CMB experiments. The extensive 
literature on this topic has been summarized in a recent article [4]. For Gaussian, statistically isotropic CMB 
sky, the Ci that correspond to covariance that maximize the multivariate Gaussian PDF of the temperature 
map, AT(q) is the Maximum Likelihood (ML) solution. Different ML estimators have been proposed and 
implemented on CMB data of small and modest size [5-10]. While it is desirable to use optimal estimators 
of Ci that obtain (or iterate toward) the ML solution for the given data, these methods usually are limited 
by the computational expense of matrix inversion that scales as N% with data size Nd [11, 12]. Various 
strategies for speeding up ML estimation have been proposed, such as, exploiting the symmetries of the scan 
strategy [13], using hierarchical decomposition [14], iterative multi-grid method [15], etc. Variants employing 
linear combinations of AT(q) such as ai m on set of rings in the sky can alleviate the computational demands 
in special cases [16, 17]. Other promising exact power estimation methods have been recently proposed [18- 
20]. 

However there also exist computationally rapid, sub-optimal estimators of Cj. Exploiting the fast spherical 
harmonic transform (~ N^ 2 ), it is possible to estimate the angular power spectrum Ci = (|a; m | 2 )/(2Z + 1) 
rapidly [21, 22]. This is commonly referred to as the pseudo-C; method [23]. (Analogous approach employing 
fast estimation of the correlation function C(q • q') have also been explored [24, 25].) It has been recently 
argued that the need for optimal estimators may have been over-emphasized since they are computationally 
prohibitive at large I . Sub-optimal estimators are computationally tractable and tend to be nearly optimal in 
the relevant high I regime. Moreover, already the data size of the current sensitive, high resolution, 'full sky' 
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CMB experiments such as WMAP have compelled the use of sub-optimal pseudo-C; related methods [26, 27]. 
On the other hand, optimal ML estimators can readily incorporate and account for various systematic effects, 
such as non-uniform sky coverage, noise correlations and beam asymmetries. 

In the years after the COBE-DMR observations [28], more sensitive measurements at higher resolution 
but with limited sky coverage were made by a number of experiments [55]. The effect of incomplete (more 
generally, non uniform) sky coverage on the sampling statistics of Ci was the dominant concern of these 
experiments such as the ground based experiment TOCO [29], DASI [30], CBI [31], ACBAR [32], and balloon 
based experiments BOOMERang [33], MAXIMA [34, 35] and Archeops [36]. Comprehensive analyzes have 
been carried out to tackle this problem. For example, the basic semi-analytic framework developed [23] was 
subsequently implemented as fast, efficient scheme for the analysis of the BOOMERang experiment [39]. 
While the non-uniform sky coverage has been addressed in the pseudo-C; method, the other effects remain 
to be incorporated. 

In this paper, we initiate a similar line of research to address a more contemporary issue that has gained 
relative importance in the post WMAP [26] (and pre-Planck) era of CMB anisotropy measurement with 
'full' sky coverage. It has been usual in CMB data analysis to assume the experimental beam response 
to be circularly symmetric around the pointing direction. However, any real beam response function has 
deviations from circular symmetry. Even the main lobe of the beam response of experiments are generically 
non-circular (non-axisymmetric) since detectors have to be placed off-axis on the focal plane. (Side lobes 
and stray light contamination add to the breakdown of this assumption). For high sensitive experiments, 
the systematic errors arising from the beam non-circularity become progressively more important. Recent 
CMB experiments such as ARCHEOPS, MAXIMA, WMAP have significantly non-circular beams. Future 
experiments like the Planck Surveyor are expected to be even more seriously affected by non-circular beams. 

Dropping the circular beam assumption leads to major complications at every stage of the data analysis 
pipeline. The extent to which the non-circularity affects the step of going from the time-stream data to sky 
map is very sensitive to the scan-strategy. The beam now has an orientation with respect to the scan path 
that can potentially vary along the path. This implies that the beam function is inherently time dependent 
and difficult to deconvolve. Even after a sky map is made, the non-circularity of the effective beam affects 
the estimation of the angular power spectrum, Ci, by coupling the I modes, typically, on scales beyond the 
inverse angular beam-width. 

Barring few exceptions (eg., [38]), the non-circularity of beam patterns in CMB experiments has been 
addressed in limited context. When it has not been totally ignored, one has measured with numerical 
simulations the biasing effect on the power spectrum of CMB anisotropics of neglecting the non-circularity 
of the beams in the data analysis chain (see e.g., MAXIMA [35, 40], Archeops [36, 37]). This approach only 
deals with the diagonal part of the matrix relating the observed power spectrum to the underlying power 
spectrum, so does not fully describe the effect of the beam complexity on the CMB statistics. An integrated 
approach to account for the systematic effect of a non-circular beam has not yet been developed. 

In this initial work we skip over the issues related to map making and focus on the CMB power spectrum 
estimation from a CMB sky map made with an effective beam that is non-circular. Mild deviations from 
circularity can be addressed by a perturbation approach [41, 42]. Besides providing an elegant analytic 
formalism, the approach has lead to rapid methods for computing the window functions for CMB experi- 
ments [43]. In this work the effect of beam non-circularity on the estimation of CMB power spectrum is 
studied analytically using this perturbation approach. 

We present a brief primer on the connection between CMB power spectrum and the experimental window 
functions in section II. The section is designed to keep the paper self-contained and also serves to set the 
notation for the rest of the paper. In section II B, we briefly review the perturbation approach for computing 
the the window functions for CMB experiments with non-circular beam [41] and also define the elliptical 
Gaussian beam and its spherical transform. The bias matrix accounting for the non-circularity of the beam 
for the pseudo-C; estimator of CMB anisotropy is derived and discussed in section III. The error-covariance 
for the unbiased estimator is derived in section III. We conclude with a discussion of the results in section IV. 
An interesting exercise of fitting the WMAP beam maps with an elliptical Gaussian beam profile is presented 
in an appendix A. Details of the steps leading to our analytical results are given in Appendix B. 



Conventionally, the CMB temperature, AT(q), is expressed as a function of angular position, q = (9,(j)), 
on the sky via the spherical harmonic decomposition, 
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WINDOW FUNCTIONS OF CMB EXPERIMENTS: A BRIEF PRIMER 
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In an idealized noise free, CMB anisotropy sky map AT(q) made with infinitely high resolution, the 
angular power spectrum is given by 

c < = 21TT £ (|a ' m|2) ' (2 ' 2) 

m— — I 

where 

a lm = [ ^y; m (q)AT(q) (2.3) 



are the spherical harmonic transforms of the temperature deviation field AT(q). We introduce the scaled 
power spectrum C/ = (1(1 + l)/2w)Ci, that measures the power per logarithmic interval of angular scale, /. 
Eliminating a lm , we may write, 



d = 1(1 + 1) 



J dn^ J d£}q 2 ( AT(qi ) AT(q 2 )) Pi (qi • q 2 ), (2.4) 



8tt 2 

where we have made use of the expansion of Legendre Polynomials 

4 1 

fl(qi-qa) = ^fl ^ ^(01)^(02)- (2-5) 

m— — I 

If we assume the isotropy of the CMB sky, (AT(qi)AT(q 2 )) should depend only on q x • q 2 . Therefore, we 
can use Legendre expansion to show that, 

00 27-1-1 

(AT(qi)AT(q 2 )) = ^ C; P ; ( qi . q 2 ). (2.6) 

All CMB anisotropy experiments measure differences in CMB temperature at different locations on the 
sky. A step of map-making is required to derive the above temperature anisotropy map at each point on the 
sky. Since this is a linear operation, the correlation function of the measured quantity for a given scanning 
or modulation strategy can always be expressed as linear sum of 'elementary' correlations of the temperature 
given in eq. (2.6). 

Typically, a CMB anisotropy experiment probes a range of angular scales characterized by a window 
function Wi(q, q'). The window depends both on the scanning strategy as well as the angular resolution 
and response of the experiment. However, it is neater to logically separate these two effects by express- 
ing the window Wj(q, q') as a sum of 'elementary' window function of the CMB anisotropy at each point 
of the map [41]. In this work, we only deal with these elementary window functions. For a given scan- 
ning/modulation strategy, our results can be readily generalized using the representation of the window 
function as sum over elementary window functions (see, e.g., [41, 43]). Although the quantitative results 
we present in this paper refer to a scan strategy where each pixel is visited by the beam only once, this is 
not a limitation of our approach. If pixels are multiply visited by the beam with different orientations, the 
correlation function still can be expressed as a sum over appropriate elementary window functions for which 
all the results we describe in this paper hold. 



A. Window function for circular beams 

Due to finite resolution of the instruments, the 'measured' temperature difference AT(q) along the direc- 
tion q in response to the CMB anisotropy signal AT(q') is given by 

AT(q) = y^B(q,q')AT(q') (2.7) 

where the experimental "Beam" response function B(q, q') describes the sensitivity of the measuring instru- 
ment at different angles around the pointing direction. There is an additional contribution from instrumental 
noise denoted by n(q) which we shall introduce later into our final results. 

The two point correlation function for a statistical isotropic CMB anisotropy signal is 

C(q, q') = (AT(q)AT(q')) = £ Q W t (q, q') , (2.8) 
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where Ci is the angular spectrum of CMB anisotropy signal and the window function 

^(qi,q 2 ) = /dnq [dSl v B{$ 1 ,®B(% i ,$)Pi(§-#), 



(2.9) 



encodes the effect of finite resolution through the beam function. 

For some experiments, the beam function may be assumed to be circularly symmetric about the pointing 
direction, i.e., B(q, q') = B(q ■ q') without significantly affecting the results of the analysis. In any case, this 
assumption allows a great simplification since the beam function can then be represented by an expansion 
in Legendre polynomials as 



B(q-q') = ^ f>Z + l)5zP;(q-q') 



Consequently, it is straightforward to derive the well known simple expression 

Wi(q,q') = SfP;(q-q'), 
for a circularly symmetric beam function. 



(2.10) 



(2.11) 



B. Window function for non-circular beams 



Pole 




FIG. 1: The figure illustrates that a beam pointed in an arbitrary direction q = (9,<f>), with an orientation given 
by the angle p(q) can always be rotated to point along z oriented with p(z) = 0. The Euler angles of this rotation 
are clearly seen to be (9,<f>,p). Consequently, the beam transforms are related through Wigner rotation matrices 
corresponding to the same rotation. 



While some experiments may have circularly symmetric beam functions, most experimental beams are non- 
circular to some extent. The effect of non-circularity of the beam has become progressively more relevant 
for experiments with higher sensitivity and angular resolution. The most general beam response function 
can be represented as 



oo 

B(z, q ) = X) L W£)Wq) 



(2.12) 



1=0 m=-l 



by a spherical harmonic expansion when pointing along z axis ("North pole" in some given astronomical 
coordinate system). In case of circularly symmetric beams, the real coefficients B\ = y 'Ait / (21 + l)6;o- 

For mild deviations, the non-circularity of the beams can be parameterized by a set of small quantities 
(3i m = bim/bio - the Beam Distortion Parameters (BDP). The smoothness of the beam response implies that 
at any multipole I, the coefficients [3i m decrease sufficiently rapidly with increasing \m\. In addition, for the 
rest of paper we assume that the beam function has reflection symmetry about two orthogonal axes on the 
(locally flat) beam plane, which ensures that the coefficients bi m {i) are real and zero for odd values of m. An 
example of a non-circular beam with such symmetries is the elliptical Gaussian beam. A brief mathematical 
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FIG. 2: The left panel shows the WMAP Ql (A side) beam map overlaid with IRAF fitted ellipses over iso-intensity 
contours. More details are in Appendix A. On the right panel, we plot the product of beam distortion parameters 
for the elliptic Gaussian fit to the WMAP-Q1 beam versus multipole corresponding to the different order of the 
perturbation expansion of a window function for a non-circular beam. Note that the effect kicks in at la ~ 1. 



description of such beams can be found later in this section. In order to verify our analytical results, we 
have used the elliptical Gaussian beam as a model of non-circular beam. However, our analytic results would 
apply to a general form of non-circular beam (as long as flu is zero or sub-dominant to fin)- 

In order to find an expression for window function in terms of the fi\ m and we follow the approach 
in [41]. The beam transforms for an arbitrary pointing direction q may be expressed as, 

i 

bim(q) = ]T b lm ,(z)D l mm ,(q,p(q)), (2.13) 

m'=— I 

where D l mm , (q, p) = D l mm , (</>, 0, p) are the Wigner-D functions given in terms of the Euler angles describing 
the rotation that carries the pointing direction q to z-axis, as illustrated in Figure 1. The third angle p(q) 
measures the angle by which the beam has rotated about the new pointing direction, when the pointing 
direction moves from z to q [56]. Inserting the spherical transform of the beam in eq. (2.13) into eq. (2.9) 
we can write the window function as 

wi(qi.qa) = &L(qi)Wq 2 ) (2.14) 

m=—l 

l l l 

= Bf Yl E AmxAm 2 ]T ^ mi ( qi ,p( qi ))^ m2 (q2,p(q2)) (2.15) 
mi— — I m2 — — I m— — l 

solely in terms of the circular component of the beam function Bi and non-circular parts encoded in the 
BDP's, Pi m . As pointed out in [41], the window function expressed in the form of eq. (2.15) has an obvious 
expansion in perturbation series in /3; m retaining only the lowest values of |mi| and \m2\- In this paper, 
we adopt this perturbation approach to evaluate the leading order correction to power spectrum estimation 
arising due to mild deviations of the beam from circular symmetry. 

For numerical evaluation it is advantageous to use the summation formula of Wigner-D to combine the 
product of the two Wigner-Z? functions in eq. (2.15) into a single one as [41] 



Wi(c Ll ,q 2 ) = ^-j E [birn>(2)Yk m „(z)D L m , m „(a- Pl , 7 ,[3 + p 2 ), (2.16) 

in' — — l m" — — l 

where 

COS7 = qi • q 2 

cota = — cos6>i cot(0i — fo) + sin#i cot#2 csc(0i — (^2) (2-17) 

COt/3 = — COS#2 COt(0i — 4>2) + COt#i Sin#2 CSc(0i — (^2). 
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TABLE I: In literature, the elliptical beams have been described by several parameters which can all be expressed in 
terms of the Gaussian widths along the semi-major (<n) and the semi-minor (172) axes of the ellipse. We have used 
these parameters at several places in the paper. 



Parameter Symbol Expression 



Eccentricity e W 1 — ^§ 

CT 2 

Non- Circularity Parameter e — 1 

CT 2 

Ellipticity e 1 — ga. 



For large values of I it is computationally expensive to evaluate the entire ml and m" sum in eq. (2.16). 
However, for a smooth, mildly non-circular beam function, restricting the summation to a few low values of 
to' and to" results in a good approximation. The leading order terms in the perturbation [41] 



Wi(qx, q 2 ) = 



47T 

21 + 1 



Mz)f 4o(7) + 26 I0 (i)6, 2 (8) {cos(2(a - Pl )) + cos(2(/3 + p 2 ))} dj, 2 ( 7 ) 
+2 Mz)] 2 [cos(2(a + /3 + p 2 - Pi))4 2 (7) + (-lj-'cos^a - - Pl - p 2 )))4 2 (^ - 7)] 



+26j (z)6j 4 (z) {cos(4(a - pi)) + cos(4(/3 + p 2 ))} 4 4 ( 7 ) 



(2.18) 



can be readily evaluated using recurrence relations similar to that of Lcgcndrc function. In the above we have 
restricted to the common situation of beam functions with reflection symmetry (0i m are real and 0i m = 
for odd m) such as the elliptic Gaussian beam described next. 

An elliptic Gaussian beam profile, pointed along the z-axis is expressed in terms of the spherical polar 
coordinates about the pointing direction as follows [41] 



B(z,q) = 



1 



27T<7l<7 2 



cxp 



2a 2 (0) 



(2.19) 



where the "beam-width" a{4>) = [ct 2 /(1 + e sin 2 0)] 1 / 2 and the "non-circularity parameter" e = (a 2 /crf — 1) 
are given in terms of 0\ and <r 2 - the Gaussian widths along the semi-major and semi-minor axis, re spectively. 
However, we characterize an elliptical beam using two different parameters: eccentricity e = \J\ — a\j a\ 
and the size parameter 0i/ 2 , the FWHM of a circular beam of equal "area" [57]. 

For elliptical Gaussian beams the spherical harmonic transform is available in the closed analytical form 



hr, 



21 + 1 (l + m)\ 



47T (/ - to)! 



(I + l/2)- m I m/2 



(Z + 1/2)V 2 



exp 



(l + l/2fal 



l- 6 - 
2 



(2.20) 



where I v {x) is the modified Bessel function [41, 53]. Note, in the above equation we have used eccentricity 
e instead of the non-circularity parameter e = e 2 /(l — e 2 ) used in [41]. (Please see Table I for the various 
definitions and characterizations of elliptical beams.) 

Fig 2 shows one of the WMAP beams as an example of a distinctly non-circular beam (see iso-contours in 
the left panel) that can be efficiently handled by the leading order term in the perturbation approach (see 
the right panel). Details of the exercise of fitting elliptical Gaussian beam profile to the WMAP beam maps 
is given in appendix A. 



III. BIAS MATRIX 

Given the observed temperature fluctuations AT(q), a naive estimator for the angular power spectrum 
based on eq. (2.2) is given by 

~ /(/ + !) 1 ^ 2 

m— — I 
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where 

a lm = y"dn 4 y^(q)AT(q)i«(q) (3.2) 

are the coefficients of the spherical harmonic transform of the CMB anisotropy map [21, 22]. The weight 
function w(q) accounts for non- uniform/incomplete sky coverage and also provides a handle to weigh the data 
'optimally'. Without the inconsequential 1(1 + 1) scaling, this naive estimator is referred to as the pseudo-C; 
in recent literature [23]. The 'pscudo' refers to fact that the estimated C; is biased. Moreover, this is a 
sub-optimal estimator of the power spectrum. This naive power spectrum estimate has to be corrected for 
observational effects such as the instrumental noise contribution, beam resolution, incomplcte/non- uniform 
sky coverage. Nevertheless, the pseudo-C; method is a computationally fast and economical approach and 
is currently a method of choice for the recent large CMB anisotropy datasets (at least for large I within the 
hybrid schemes [4]). 

Faced with the computational challenges of large data sets, an approach that has been adopted is to 
compute the pseudo-C; 's from the CMB observations and then correct for the observational effects. The 
true C; spectrum is linearly related 

(Ci) = ^2A w C v (3.3) 

v 

to the pseudo-C; through a bias matrix Aw . Similar bias matrices arising due to the effect of non-uniform sky 
coverage, instrumental noise have been studied [23, 39]. In this paper, we compute the Aw for non-circular 
beam and give explicit analytical results for the leading order terms for non rotating beams. 
The pseudo-C; estimator in eq. (3.1) can be expressed as 

C; = J dfi 4l J dfi^ ™(qiMq 2 ) AT( qi )AT(q 2 )P ; ( qi • q 2 ). (3.4) 

The ensemble expectation value of the pseudo-C; power spectrum estimator is 

^ ^Tif '^^'^ / dn * I d ^ ,B (^^) B (^ q')fl'(q-q')- (3-5) 

Recalling the definition of a window function in cq. (2.9), the most general form of the bias matrix 

jmr^Y) J dn ^J daWqiMq 2 )^qi-q2)^(qi,q 2 ). (3.6) 



21' + 1 1(1 + 1) 

"' ~ 16^ 2 l'(l' + 1) 

Using the expression for the window function for a non circular beam in eq. (2.15) the bias matrix can be 
written as 



Bf, (2V + 1) 1(1 + 1) X 

'"- 47 r(2/ + l)^' + l) X i;A 



^'m' /^q^(q)^w(q,p(q))^(q 



2 



(3.7) 



The above expressions in eq. (3.6) and eq. (3.7) are valid for a completely general non-circular beam with an 
arbitrary orientation at each point. The scan pattern of the CMB experiment and relative orientation of the 
beam along it is encoded in the function p(q). The weight w(q) can account for non-uniform sky coverage. 
Analytical progress can be made when w(q) = w(9) and p(q) = p(9) are fixed along a given declination, 
but we do not discuss further it here. When the beam transform, weight function and the scan pattern 
are specified, the bias matrix can be evaluated numerically using eq. (3.7). However, for mild deviations 
from circularity, the above expression also points to a perturbation expansion in the small beam distortion 
parameters, /3; m . 

For obtaining fully analytical results, we set the weight function w(q) = 1, corresponding to a full, uniform 
sky coverage and also limit attention to scans with 'non-rotating' beams where p(q) = 0. This is presented 
in the next subsections 
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A. Circular Symmetric Beam 

We first consider eq. (3.7) for the simpler and well studied case of a circular beam. For clarity of presen- 
tation, we limit our discussion full, uniform sky coverage (u>(q) = 1). Results for non-uniform coverage with 
a circular beam are available in the literature [4, 23, 39]. 

Using the expression for the window function for circular beam eq. (2.11) into the expression for the bias 
in eq. (3.6) we recover 

Aw = Bf5 w => (Ci) = BfC t . (3.8) 

For a full sky measurement with a circular beam, the bias matrix is diagonal implying that there is no 
mixing of power between different multipolcs. The true expectation value of the power spectrum can be 
obtained by dividing the pseudo-C; estimator by the isotropic beam transform Bf. 

Next we account for the noise contribution and recover the well known result for a full sky observation. 
The pixel noise n(q) adds to the observed temperature, so that the resultant observed temperature 



AT(q) = AT(q) + n(q) (3.9) 

and we can readily obtain 

{C[) = (C t ) + = BfC t + Cf, (3.10) 

where Cf 1 is the angular power spectrum of the noise n(q) is a well determined quantity. The unbiased 
estimator for Ci obtained is 

cf B = s r 2 (5 - c?) . (3.ii) 



B. Non-circular Beam 



We obtain analytic results for the bias matrix for a full sky observation (u>(q) = 1) with a non-circular 
beam that 'does not rotate'. The phrase "non-rotating" means that the orientation of the non-circular beam 
does not rotate about its axis (the pointing direction) while the pointing direction scans the sky implying 
that 

p(q) - 0. (3.12) 

For non-rotating beam, the calculation of the bias is completely analytically tractable. The integral in the 
expression for the bias in eq. (3.7) is given by 

J i£(q) Ojdfy = y/(2l + l)irl% m ,5 mn , (3.13) 

where 

d L mO (0)d L mml (6)dco S 9, (3.14) 

and d l mm , (9) are Wigner-c? functions related to Wigner-D functions 

D l mm ,(q,p) = e-" n U l mm ,(8)e-" n 'P. (3.15) 

The analytic simplicity arises from the fact that for p(q) = 0, the Wigner-D function reduces to spherical 
harmonic function 



= \l^-D l m0 (q,0). (3.16) 



In deriving the above have used the orthogonality of the phases e n ^d(f) = 2ttS„ 
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Substituting the expression for the integral cq. (3.13) into the expression for the bias in eq. (3.7), we obtain 

2 



m——L 



V 



^1 Pl'm' I mm' 



m'=—V 



(3.17) 



where L = min{Z,Z'} is the smaller between Z and I'. 

Further analytical progress is possible for smooth beam with mild deviations from circular symmetry 
through a perturbation in terms of the small beam distortion parameters, [3i m . We calculate the exact 
analytic expression for the leading order effect. Assuming a beam with reflection symmetry where (3i m are 
zero for odd m, the leading order effect comes at the second order, namely, flnftm ( see ec L- (2-18)). Neglecting, 
film for |m| > 2, we obtain 



2 (2V + 1\ 1(1 + 1) A r lV , if 11' 



(3.18) 



Next we obtain analytical expression for the two integrals, and I^ 2 + /^_ 2 . The first one can be found 
in standard texts (e.g. [44]) given as 



jll 1 = f 1 

'mO — J 



d l mQ (e)d l mQ (6)d cos e = 



21 + 1 



$11' ■ 



(3.19) 



For m = 0, writing d l 00 (9) and d l Q2 (0) in terms of P;(cos#) and its first derivative P/(cos#) we have shown in 
Appcndix-B that for odd values of Z + Z', I$ 2 + Iq_ 2 = 0. For even values of I + I', 



jll' , jll' _ 1 

J 02 ' J 0-2 — S 



S/k if Z < Z' 

if I > I', 

{ -(4Z//t)(Z- 1)/(2Z + 1) ii 1 = 1' 



(3.20) 



where k = - 1)1' (I' + 1)(Z' + 2). 

To evaluate I^ 2 + I l J n _ 2 for non-zero m we expand d l m±2 {9) in terms of d l m0 (9) using a recurrence relation 
of the Wigner-£> functions (where I" takes integer values between I' — 2 to I' + 2) . The details are given in 
Appendix-B. We obtain that I l Jn 2 +Im- 2 = for odd l+l'. For even values of l+l', if L = min{Z, I'} > \m\ > 0, 



'(4/*)(N+i)>/SiKS in<r 



jll I T' 



ra2 T -i m _2 



(3.21) 



_ (4//s)[|m| - (Z 2 + Z + 1)/(2Z + 1)] if Z = Z'. 

The bias matrix including the leading order beam distortion (for non-rotating, reflection symmetric beams) 
can be summarized as 
• For odd values of l + l', 



Aw 



0. 



(3.22) 



• For even values of I + I' , 



( (Ft, R, ,\ 2 ( .. 8'('+l)(2'' + l) ) [o 1 (l+my.(l'-m)\ ( , n 2 

\Dl'PV2) ^'2(;'+l)2(V-l)(;'+2) ) |_ Z + 2^m=l (l- m y.(l'+m)\y m + 1 ) 



(R,,Ft„ .) 2 ( ... 8 ' ( ' +1)(2 ' , + 1) IV'' C-m)!q'+m)! r v 

(&l'Pl>2) [l'i(i' + 1 )i(i^i ){ i' +2) J [l^ m =l (l+ m )\(l'-m)\\ m ~ l > 



21+1 



(i 2 +i + l)-(2i+l)m 

V(;-i)'('+i)('+2) 



if z < /' 
if Z > /' 

if I = I'. 



(3.23) 
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FIG. 3: The log of normalized bias matrix An* /(BiBii) is plotted for an elliptical beam of eccentricity e = 0.6 and 
mean beam-width a = 0.074. The normalization is carried out so that the effect of non-circularity on the bias matrix 
can be easily compared to that for circular beams. Beam rotation and cut-sky effects have not been considered in 
this figure. One notices that the off-diagonal elements of the bias matrix take significant values for la > 1 . 

The non-zero off-diagonal terms in the bias matrix Aw imply that the non-circular beam mixes the 
contribution of different multipoles from the actual power spectrum in the observed power spectrum. Off- 
diagonal elements in Aw that arise from non- uniform/incomplete sky coverage have been studied earlier and 
are routinely accounted for in CMB experiments. Non-circular beam is yet another source of off-diagonal 
terms in the bias matrix and should be similarly taken into account. In general, CMB experiments have 
both non-circular beams and non- uniform/incomplete sky coverage that could lead to interesting features in 
A w - 

Although the analytical result is limited to mildly non-circular and non-rotating beam functions, it does 
bring to light certain generic features of the effect of non-circular beam functions. To be specific, we compute 
the elements Aw for non-rotating elliptic Gaussian beams (see appendix A) . The non-circularity of these 
beams is characterized by their eccentricity e = — o\I<j\, where <j\ and a-i are the la beam-widths along 
major and minor axes of the beam (see table I). Many experiments have characterized their beams in terms 
of an elliptic Gaussian fit (e.g. ,[40, 42, 43]). A convenient advantage of elliptical beams is that the beam 
transform bi m (and obviously, the beam distortion parameters, (3i m ) can be expressed in a closed analytical 
form. The results expressed in terms of la are broadly independent of the average beam-size [41]. 

Fig. 3 shows a density plot of the normalized bias matrix Aw / '(BiBy) for a non-rotating elliptical beam. 
The plot illustrates the importance of off-diagonal terms that arise due to the non-circular beam relative to 
the diagonal terms. The absence of coupling between multipoles separated by odd integers is evident. Also 
evident is the fall off as one moves away from the diagonal. The left panel of Fig. 4 shows that the off-diagonal 
elements of Aw are important at la ~ 1. The results are qualitatively independent of the average beam size 
a. The right panel Fig. 4 shows the strong dependence of the dominant off-diagonal element Au+2 on the 
eccentricity of the beam. 

The analytical results and numerical computations using eq. (3.6) were compared. The numerical and 
analytical results match perfectly as shown in Figure 4. Numerical computation involves the pixelized sky and 
the algorithm must ensure that this does not introduce spurious effects. We verify that Aw has numerically 
negligible off-diagonal elements when the beam is circularly symmetric. The numerical computation for 
non-circular beam are verified to be robust to the pixelization of the sky. 

Next we illustrate effect of beam-rotation and non-uniform sky coverage for a hypothetical experiment 
where Aw have been computed numerically. The left panel of Fig. 5 shows (in log scale) the normalized bias 
matrix arising from a 2.5° circular beam including a non trivial io(q) in the form of a smoothed version of 
the galactic mask Kp2 of WMAP [26, 45]. The right panel of the figure shows the extra effect that a rotating 
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non-circular beam would introduce. We assume a simple 'toy' beam rotation along an equal declination scan 
strategy, where the beam continuously 'rotates' by 2tt for every complete pass at a given declination which 
implies the simple form 



The elements here have been computed numerically using eq. (3.6) retaining the leading order terms in the 
perturbation expansion of Wi in eq. (2.18). The off-diagonal effects at low / are dominated by the cut sky 
effect. The off-diagonal element la > 1 arise solely due to non-circular beam. The numerical computation 
illustrates the potentially large corrections that can arise due to non-circular beam that 'rotate' on the sky. 
The numerical computations in this work pave the way for introducing realistic scan-pattern, beam-rotation 
and non-uniform sky coverage in a future extension to our work. 



FIG. 4: Elements of the bias matrix Aip are plotted in this figure as a function of multipole (I). The bias matrix 
relates the observed Cis to their true values. When non-circular beams are used in CMB experiments, the bias matrix 
can be shown to be non-diagonal, thus implying mixing of power between multipoles. On the left panel, we plot A u i 
for l' — I = 2, 4, 6. It is evident that the effect decreases as we move away from the diagonal and that it kicks in at 
la ~ 1, for a beam of eccentricity e = 0.8. For the figure in the right panel, we plot Au+2 for several beams of the 
same size but different eccentricities. Clearly, the effect also depends strongly on the non-circularity of the beam. 

We summarize the following features of the bias matrix : 

1. There is no coupling between (Ci) and for odd values of I + I', 

2. Coupling decreases as |Z — l'\ increases, 

3. Coupling increases with eccentricity for fixed beam size, and 

4. Size of the beam determines the multipole I value for which coupling will be maximum (la ~ 1). 

Figure 6 roughly indicates the level and nature of the effect of neglecting the non-circularity of the beam 
on CMB power estimation (for the conservative case of non-rotating beams). Consider the power spectrum 
Ci = y~)[, AwCi> measured using a non-circular, elliptical Gaussian beam of a given eccentricity, e and average 
beam- width, a. We compare the power spectrum obtained by deconvolving Ci with a circular, Gaussian beam 
of the beam-width, a with the true Ci. The lower panel shows that the error can be significant for multipole 
values beyond the inverse beam- width even for modestly non-circular comparable to the WMAP beam maps 
(Q-band) discussed in the Appendix A. 

Finally, we construct the unbiased estimator for the angular power spectrum. Invoking steps similar to 
the case of circular beams to account for the instrumental noise, we obtain 



p(q) = p(0,<£)=0. 



(3.24) 




= 0.2 Analytic 
= 0.2 Numeric 
= 0.4 Analytic 
= 0.4 Numeric 
= 0.S Analytic 
= 0.B Numeric 
= 0.8 Analytic 
= 0.8 Numeric 




(3.25) 



The unbiased estimator for the angular power spectrum is 




(3.26) 
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FIG. 5: The normalized bias matrix elements (on log scale, \og[Aur /(BiBy)]) of a hypothetical experiment with a 
scan pattern (eq. 3.24) corresponding to a rotating, non circular beam (e = 0.6) and non uniform sky coverage are 
studied. The left panel shows the effect of non-uniform coverage alone (circular beam approximation). The right 
panel isolates the additional effect that arises due to the non-circularity of the beam and its rotation. We note that 
significant off-diagonal elements arise at la > 1 from the non-circular beam comparable to that from the non-uniform 
coverage. The non-uniform coverage corresponds to a smoothed WMAP Kp2 galactic mask (smoothed from resolution 
of N 3 ide = 512 to 64). We use a sufficiently high resolution beam with a = 0.018 (#1/2 = 2.5°) to ensure that the 
effects due to the galactic mask and the non-circular beam appear in distinct regions of the multipole space. 
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FIG. 6: The effect of non-circular beam is studied for CMB power spectrum estimation by a CMB experiment with 
a WMAP-like non-circular beam. For illustration, we consider the best fit (Power law) model to be the (fiducial) 
true Ci of the Universe shown as the solid line in the upper panel. Let Ci be the power spectrum measured by using 
a elliptical, Gaussian beam with eccentricity, e = 0.6, and a — 0.0016. The dashed line shows the Ci that would be 
inferred by deconvolving Ci with a circular beam assumption with beam-width, a. The lower panel plots the relative 
error in the power spectrum recovered with a circular beam assumption for a measurements made with a non-circular 
beam with e = 0.4 to 0.8. 
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IV. ERROR-COVARIANCE MATRIX 

The statistical error-covariance of the estimated angular power spectrum is defined as 

Cov(C,,C,0 = {{Ci-{Ci)){C v -{C v ))). (4.1) 
In an idealized, noise free, CMB experiment with infinite angular resolution uniformly covering the full sky 

& = / J ^q 2 AT(q 1 )AT(q 2 )P i (q 1 • q 2 ). (4.2) 

Using the property of Gaussian random fields that, 

(AT(q 1 )AT(q 2 )AT(q' 1 )AT(q' 2 )) = (AT(q 1 )AT(q 2 ))(AT(q' 1 )AT(q' 2 )) + (4.3) 

(AT(q 1 )AT(q / 1 ))(AT(q 2 )AT(q 2 )) + (AT(q 1 )AT(q 2 ))(AT(q / 1 )AT(q 2 )) 

and eq. (2.6), we recover the well known result for full sky CMB maps 

CoviCuCv) = ^^(CtfSu, = ^jCfS u ,, (4.4) 

corresponding to Ci being a sum of the squares of 21 + 1 Gaussian variates, i.e. X21+1 distribution. The 
measured power spectrum at each multipole is independent (for full sky CMB maps). The variance of the 
power spectrum estimator is not zero even in the ideal case. Consequently, the measurement angular power 
spectrum from the one available CMB sky map is inherently limited by an inevitable error the Cosmic 
Variance [58]. 

A. Circular Beam 

For measurements made with a circular beam, the temperature is a linear transform of the actual temper- 
ature (see eq. (2.7)). So, it also represents a Gaussian random field. Hence, eq. (4.3) remains valid even for 
observed temperature fluctuations. Moreover, the window function takes a simple form given in eq. (2.11). 
Consequently, eq. (2.6) gets modified to 

* 2Z + 1 



(AT(qi)AT(q 2 )) = £ Bf d fl(qi • q 2 ). (4.5) 



The covariance matrix 



remains diagonal for circular beams, i.e., the measured power spectrum at each multipole is independent of 
the power measured in the other multipoles. The second equality follows from eq. (3.8). 
Including the instrumental noise spectrum in the measured power spectrum Cf , we obtain 

Covert) = ^((C^+C^) 2 , (4.7) 

where we assume that the noise spectrum Cf is known much better and, in particular, does not suffer from 
cosmic variance. For the unbiased estimator given by eq. (3.11), the well known covariance matrix 

Cov(Cr,c7 B ) = Sr 4 Cov(C~',ClO = ^(Ct+B^Cn 2 (4.8) 
is readily obtained from the linear transformation between C[ and C\ B [48, 49] . 
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B. Non-circular Beam 



As expected, the covariance for the non-circular beam is considerably more complicated. We start with 
the general form of the two point correlation function. Using eq. (3.4), the general form of the covariance 
matrix is 



Cav(Ci,Ci') = 



(4^) 4 



E 

h,h=0 



(2fi + l)(2f 2 + l) 

Wi + i)(i 2 + i) 



C h Ci 2 / d 4 nw(qi)w(q 2 )w{c{ 1 )w(q / 2 ) x 



fl(qi • q2)^'(q'i ■ q' 2 ) [^(qi,qi)^ 2 (q2,q 2 ) + W h (qi, <QW h (q 2 , qi)], 



(4.9) 



where for brevity we denote d 4 fl = dil^ dQ,q 2 dfl^ dCl^ . 

Noting the interchangeability of the dummy variables q' x and q' 2 , we combine the two terms in the above 
equation to obtain 



Cav{Ci,C v ) = 2 



(4^) 4 



h,h=0 



d 4 fMqiMq 2 MqiMq 2 )^(qi • q 2 )P/'(qi • qiW^qi.qiW^qa.qi). (4.10) 

We expand the Legendre Polynomials in terms of spherical harmonics (eq. (2.5)) and use the expression for 
the window function in eq. (2.15) to obtain 

Cav{C u Cv) - - ^ - - C h C h B h B l2 ^ ^ ^ ^ 

(i,( 2 =0 m= — ( m' = — (' mi=— ii m 2 =— (2 

E /WAk / dn« I «'(qi)^(qi)^ iroi (qi,p(qi)) / ^^(q;)^^,^)^^^,^)) 

E A^A^ I dOq 2 «;(q 2 )y ;m (q 2 )^; m ,(q 2 ,p(q 2 )) J dfi^ w^)Y{! m , (q 2 )£>£ 2ro ,, (q' 2 , p(q' 2 )) 

m 2 ,m 2 ' = — Z2 

(4.11) 

as the general expression for error covariance for angular power spectrum for non-circular beams. Note that 
even for full, uniform sky observations, w(q) = 1, the error covariance matrix is no longer diagonal. 

To make further progress analytically, we restrict to the case of uniform, full sky coverage (w(q) = 1) with 
no beam rotation (p(q) = 0). Using the integration of eq. (3.13) and after a considerable algebra we may 
write the expression for covariance as 



Cov(C,,a,, = "' ( ' +1 '"' + " £ 



m— — L 



I I 1 2 

l\ — \m\ ' m^ = — Zi m"=— /i 



(4.12) 



where L — min{Z,Z'} is the smaller between I and I'. The integrals 7^' m / are defined in §111 and the 
analytical expressions for m! — 0, ±2 are given. It is straightforward to verify that the above equation 
correctly reproduces the expression for the error-covariance in the circular beam case given by eq. (4.6). 

For evaluation of the covariance matrix, we note that though the summation over l\ runs from to oo, the 
contributions are significant only around I <~ 1/a and the summation can be truncated suitably. Further, for 
most beams we can confine to the leading order approximation as in eq. (2.18), by neglecting all the /3; m 's 
for m > 4. For mild deviations from circular beams, the observed power spectrum at different multipolcs 
are weakly correlated ( ~ PnPi'2)- The error-covariance matrix can be diagonalized to find the independent 
linear combinations of estimators (eigenvectors), and the variances of theses independent estimators are 
given by the corresponding eigenvalues. These eigenvalues are necessarily larger that the cosmic variance 
corresponding to a circular beam. 

The inclusion of instrumental noise is similar to what was done in the circular beam case. The covariance 



cov(c;,c;o = cov(g,c,o + 



2S W 
21 + 1 



2lfi l )C? + (C/ v ) 



N\2 



(4.13) 
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FIG. 7: Log of the normalized covariance matrix Cov(C;, CV )/(-£?;£?;< ) 2 [in the units of (fiK) 4 ] is plotted for an 
elliptical beam of eccentricity e = 0.6 and mean beam-width a — 0.074. Due to the non-circularity of the beam, the 
error in CMB angular power spectrum estimate at different multipoles are no longer independent. We notice that 
the off-diagonal elements of the error covariance matrix are pronounced for la > 1. 



clearly reproduces the result in eq. (4.6) in the limit of a circular beam. Figure 7 shows a density plot of the 
elements of the covariance matrix for a non-circular (elliptical) beam with no rotation. In contrast to the 
case for incomplete (cut) sky case, where the effects are at small I (see [4]), the non-circular beam affects the 
large multipoles region (la > 1). The pseudo-C; approach is close to optimal for large I hence it may be more 
important to account for non-circular beams effects than the cut-sky, since it is possible to use maximum 
likelihood estimator for small /. 

The error-covariance matrix for the unbiased estimator eq. (3.26) for non-circular beams is given by 

Cov(cTV7 B ) = EEATM^Cov^,^) 

h h 

h h 

where the matrix aw = B,7 i™ , being very close to identity, demonstrates that the beam-modified cosmic 
variance part of the covariance of unbiased estimator weakly depends on Bi 's, whereas the noise part depends 
on them significantly. 



, (4.14) 



V. DISCUSSION AND CONCLUSION 



We present an analytic framework for addressing the effect of non-circular experimental beam function 
in the estimation of the angular power spectrum Ci of CMB anisotropy. Non-circular beam effects can be 
modeled into the covariance functions in approaches related to maximum likelihood estimation [9, 10] and 
can also be included in the Harmonic ring [16] and ring-torus estimators [17]. The latter is promising since it 
reduces the computational costs from iV 3 to iV 2 . However, all these methods are computationally prohibitive 
for high resolution maps and, at present, the computationally economical approach of using a pseudo-C/ 
estimator appears to be a viable option for extracting the power spectrum at high multipoles [4]. The 
pseudo-C/ estimates have to be corrected for the systematic biases. While considerable attention has been 
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devoted to the effects of incomplete/non- uniform sky coverage, no comprehensive or systematic approach is 
available for non-circular beam. The high sensitivity 'full' (large) sky observation from space (long duration 
balloon) missions have alleviated the effect of incomplete sky coverage and other systematic effects such as 
the one we consider here have gained more significance. Non-uniform coverage, in particular, the galactic 
masks affect only CMB power estimation at the low multipoles. Recently proposed hybrid scheme promotes 
a strategy where the power spectrum at low multipoles is estimated using optimal Maximum Likelihood 
methods and pseudo-C; are used for large multipoles. 

We have shown that non-circular beam is an effect that dominates at large / comparable to the inverse 
beam width. For high resolution experiment, the optimal maximum likelihood methods which can account 
for non-circular beam functions are computationally prohibitive. In implementing pseudo-C; estimation, 
the non-circular beam effect could dominate over the effects of more well studied effect of non-uniform sky 
coverage. Our work provides a convenient approach for estimating the magnitude of this effect in terms of 
the leading order deviations from a circular beam. The perturbation approach is very efficient. For most 
CMB experiments the leading few orders capture most of the effect of beam non-circularity. The perturbation 
approach has allowed the development of computationally rapid method of computing window functions [41]. 
Our work may similarly yield computationally rapid methods correcting for beam non-circularity. 

The quantitative estimates of the off-diagonal matrix elements of the bias and error-covariance for 'non- 
rotating' beam graphically illustrate the general features that can be gleaned from our analytic results. They 
show that the beam non-circularity affects the C/ estimation on multipoles larger than the inverse beam 
width. A strong dependence on the eccentricity of the beam is also seen. We caution against interpreting 
these results as a measure of the non-circular beam effects for any real CMB experiment. The analytical 
results are limited to non-rotating beams and uniform sky coverage. Numerical results do not include scan 
pattern of any known experiment. Numerical calculations of the bias matrix for a 'toy' scanning strategy 
where the beam rotates on the sky indicates the possibility of significant corrections. The bias due to non- 
uniform sky coverage can have interesting coupling to the bias from beam non-circularity. On the other 
hand, it has also been shown that effects of non-circular beams can be diluted if the scan pattern is such 
that each point in the sky is revisited by the beam with a different orientation at different time [40]. The 
numerical implementation of our method can readily accommodate the case when pixels are revisited by 
the beam with different orientations. Evaluating the realistic bias and error-covariance for a specific CMB 
experiment with non-circular beams would require numerical evaluation of the general expressions for Au> 
in eqs. (3.7) using real scan strategy and account for inhomogeneous noise and sky coverage. We defer such 
an exercise to future work. 

It is worthwhile to note in passing that that the angular power Ci contains all the information of Gaussian 
CMB anisotropy only under the assumption of statistical isotropy. Gaussian CMB anisotropy map measured 
with a non-circular beam corresponds to an underlying correlation function that violates statistical isotropy. 
In this case, the extra information present may be measurable using, for example, the bipolar power spec- 
trum [50]. Even when the beam is circular the scanning pattern itself is expected to cause a breakdown 
of statistical isotropy of the measured CMB anisotropy [39] . For a non-circular beam, this effect could be 
much more pronounced and, perhaps, presents an interesting avenue of future study. 

In addition to temperature fluctuations, the CMB photons coming from different directions have a random, 
linear polarization. The polarization of CMB can be decomposed into E part with even parity and B part 
with odd parity. Besides the angular spectrum Cf T , the CMB polarization provides three additional spectra, 
Cf E , C EE and C BB which are invariant under parity transformations. The level of polarization of the CMB 
being about a tenth of the temperature fluctuation, it is only very recently that the angular power spectrum 
of CMB polarization field has been detected. The Degree Angular Scale Interferometer (DASI) has measured 
the CMB polarization spectrum over limited band of angular scales in late 2002 [51]. The WMAP mission 
has also detected CMB polarization [52]. WMAP is expected to release the CMB polarization maps very 
soon. Correcting for the systematic effects of a non-circular beam for the polarization spectra is expected 
to become important soon. Our work is based on the perturbation approach of [41] which has been already 
been extended to the case of CMB polarization [42]. Extending this work to the case CMB polarization is 
another line of activity we plan to undertake in the near future. 

In summary, we have presented a perturbation framework to compute the effect of non-circular beam 
function on the estimation of power spectrum of CMB anisotropy. Wc not only present the most general 
expression including non-uniform sky coverage as well as a non-circular beam that can be numerically 
evaluated but also provide elegant analytic results in interesting limits. In this work, we have skipped over 
the effect of non-circular beam functions on map-making step. In simple scanning strategies, our results 
may be readily applied in this context. As CMB experiments strive to measure the angular power spectrum 
with increasing accuracy and resolution, the work provides a stepping stone to address a rather complicated 
systematic effect of non-circular beam functions. 
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APPENDIX A: ELLIPTICAL GAUSSIAN FIT TO THE WMAP BEAM MAPS 



We briefly describe an exercise in characterizing non-circular beams in CMB experiments using the beam 
maps of the WMAP mission. We analyzed the WMAP raw beam images in the Ql, VI and Wl [45, 46] 
bands using two different standard software packages. We use the elliptical Gaussian fit allowed by the 
well known radio-astronomy software, AIPS and a more elaborate ellipse fitting routine available within the 
standard astronomical image/data processing software IRAF. The ELLIPSE task in the STSDAS package of 
IRAF, which uses the widely known ellipse fitting routines by Jedrzejewski [47], allows independent elliptical 
fits to the isophotes. This significant greater degree of freedom in fitting to the non-circular beam allows 
us to assess whether a simple elliptical Gaussian fit is sufficient. The three bands see Jupiter in the two 
horns (labeled A and B) as a point source. The fitting routine fits ellipses along iso-intensity contours of the 
beam image, parameterized by position angle (PA), ellipticity (e) and position of the center. Each of these 
parameters can be independently varied. The distance betw een successiv e ellipses can also be independently 
varied. The eccentricity e is related to ellipticity e as e = yjl — (1 — e) 2 (Please see Tabic I). 




FIG. 8: The beam profile is characterized by intensity along the semi major axis (SMA). The beam in Q-band for 
WMAP experiment was analyzed using IRAF and fitted to both circular and elliptical profiles. We have plotted the 
best fit circular profile (solid error bars) and overlaid the profile recovered by inverting the WMAP beam transforms, 
available at LAMBDA website (solid line). Two analytical models for circular beam profile g(0) and h(9) are also 
considered, and the best fit profiles are overlaid. We find that these models are consistent with the IRAF and WMAP 
data. We have also plotted the best fit elliptical profile along SMA (broken error bar). Notice that the error bars in 
this case are much smaller than those for circular profile, implying a better agreement with the data. 



We fit the the beams in two different ways: (a) by holding the ellipticity constant to e = 0.05 and freely 
varying the position angle and center and (b) fixing the center to be the pixel with the highest intensity 
(normalized to 1.0 at the central pixel) and varying ellipticity and position angle. In the first case, we get 
the closest approximation to circular beam profiles as used in WMAP data analysis. This beam has no 
azimuthal (0) dependence. In the latter case, we get the elliptical profile of the WMAP beam which depends 
on both the polar (9) and azimuthal (</>) distance from the pointing direction. Notice that in this case it is 
sufficient to provide the intensity along a particular direction (usually, the semi-major axis or <f> = 0) and 
the ellipticity e. 

Even a visual inspection reveals that the Ql beam map plotted in Fig 2 is non-circular and the iso- 
intensity contours distinctly elliptical. Thus it comes as no surprise that the error bars as shown in Figure 8 
for circularized beam are larger than those for the elliptical profile. As a consistency check, we take the 
WMAP Ql beam transfer function Bi from WMAP first year data archived at publicly available LAMBDA 
site [45] and 'recover' the circular beam profile B(0) using eq. (2.10). 

From Figure 8, it is clear that this 'recovered' beam profile is in good agreement with that obtained by 
IRAF. This allows us to make some statements about the profile fitting in CMB experiments, in the context 
of WMAP beams. The beam profile B{9) has been modeled as a Gaussian times a sum of even order Hermite 
polynomials (H 2n ) by the WMAP team [46]. To compare, we have also modeled the beam profile with a 
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FIG. 9: The beams for WMAP experiment in three bands Q, V and W for both the horns (A and B) were fitted to 
elliptical profiles using IRAF. The plot above shows the fitted eccentricity and position angle along the semi major 
axis (SMA). The presence of side-bands in the smaller beams (W band) makes it difficult for IRAF to model them 
sufficiently well. However, in the Ql band, such sub-structures in the beam are not present thus allowing the IRAF 
ellipse fitting routine to fit reasonably good ellipses which have consistent eccentricities (e ~ 0.65) and position angles 
all along the SMA. The VI beam is smaller in extent than the Ql beam and its eccentricity was determined to be 
e ~ 0.46 using IRAF. The highest resolution beam in in Wl band, whose eccentricity was determined to be e ~ 0.40. 



function h(9) given by 

h(9) = cxp ^~ a9 2 ^j (h + h 2 H 2 (8) + h 4 H 4 (6)) , (Al) 

where a,h ,h 2 and h 4 are unknown parameters to be fixed by least squared method. We found that this 
model fits the data very well with a reduced x 2 of about 0.7. However, on closer analysis, it is found that 
the chief role of the Hermite polynomials is to add a constant baseline over and above the Gaussian. To test 
this hypotheses, we choose another form of the fitting function g(9) given by 

9(0) = ,9o + 9i cxp (-^ g 2 8 2 ^j , (A2) 

where go, Si and g 2 are parameters of the model. It is very interesting to note that this model also fits the 
data very well with a reduced % 2 of about 0.8 for the best fitted parameters. In all fairness, g(ff) serves 
as a simpler model for the beam profile. We cannot point to the precise origin for the baseline. However, 
such 'skirts' in beam responses are not uncommon in radio-astronomy. At this point, our observation should 
perhaps merit a curious aside, if not as an alternative approach to beam modeling. Our best-fit models g(9) 
and h(9), along with the IRAF fitted data points to the WMAP Ql (A) beam is shown in Figure 8. 

As shown earlier in this paper, the effects of non-circularity of the CMB experimental beams show up in 
the power-spectral density estimates through the off-diagonal elements of the bias matrix Aw ■ As shown in 
eq. (3.7), these in turn can be expressed in terms of the leading components of the harmonic transform of 
the beam. In general the harmonic decomposition of a non-circular beam may have to be done numerically. 
But for the particular case of an elliptical Gaussian beam, a closed form expression given by eq. (2.20) serves 
as a useful test-bed for us. Thus another motivation for fitting ellipses to WMAP beams using IRAF was 
to get a handle on the eccentricity of these beams so as to find the harmonic transform components of an 
elliptical Gaussian beam of similar eccentricity. This allows us to give more realistic estimates of the effect 
of non-circularity of the beam on Ci estimates. 
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TABLE II: The result of ellipse fitting using IRAF on the Ql, VI and Wl beams of the WMAP experiments. The 
frequency quoted is the 'effective' frequency of the corresponding band from Page et. al. [46]. The presence of 
sub-structures in the Wl band makes it difficult to fit elliptical contours to the beam. 



Beam Frequency Eccentricity Position Angle 



(GHz) 




(degree) 


Ql (A) 40.9 


0.65 


+80 


Ql (B) 40.9 


0.67 


-80 


VI (A) 60.3 


0.48 


+60 


VI (B) 60.3 


0.45 


-60 


Wl (A) 93.5 


0.40 


a 



"Not well determined by the IRAF ellipse fitting routine. 



It is interesting to note how the fitted eccentricities vary as a function of the distance along the semi-major 
axis of the fitted ellipses for various beams. The smaller beams (VI and Wl) have sufficient sub-structure in 
the form of side-lobes which throws the ellipse fitting routine off course. However, where the sub-structure 
is less pronounced, we find that the eccentricities of the fitted ellipses takes a constant value. Toward the 
center of the ellipses, there are far too few pixels to average over, which in turn manifests as large error 
bars in the eccentricities and position angles of the ellipses. In Figure 9, we notice that the Ql beam has 
a very elliptical profile with eccentricity e > 0.65 and position angle of about 75°. We also fitted the Ql 
(A) beam to an elliptical Gaussian model using radio astronomy standard data analysis software AIPS and 
got consistent numbers for the eccentricity. However the IRAF modeling gives us more freedom to vary the 
eccentricity and position angle as we move away from the center of the ellipse and the result is that the 
beam is modeled more accurately. 



APPENDIX B: DETAILS OF ANALYTIC DERIVATIONS 



In the appendix we provide the details of the analytical steps involved in deriving some of the expressions 
used in the main text. This is designed to keep the paper self contained and easy to extend. 

First, we outline the steps involved in evaluating the integral Iq 2 + ^0-2 = 2 d l 00 (9)d l 02 (9)d cos 9 to 
obtain the result in eq. (3.20). 

Using the expressions [41, 44] for d l aQ and d 02 in terms of Legendre Polynomials and its derivatives, 

/ d l Q0 (9)d l m (9)d cos 9 = - l JLl til f P l ( x )P l ,(x)dx + - [ xP l {x)P[,{x)dx, (Bl) 

7-i K 7-i K 7-i 

where, k = \JJl — 1)1(1 + 1)(Z + 2). The first integral is simply the orthogonality of Legendre polynomials 

28u> 



Pl(x)Pl,(x)dx , )j + y 



Further, we can show that for odd values of I + I', 

J xPi(x)P[,{x)dx = 0, 

and for even values of I + I', 

J xPi{x)Pl,(x)dx = < 
Assembling all these we can derive eq. (3.20). 



(B2) 



(B3) 



2 if I < V 
if I > V 
22/(22 + 1) if 1 = 1'. 



(B4) 
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Next we evaluate the more general integral I% 2 + ^m-2 = f\ d l m o(9) [<Ci2^) + ^m-2(^)] ^ cos $ to obtain 
the expression in cq. (3.21). The first step is to express dj 7l±2 (^) m terms of d m0 (6). Using the recurrence 
relations for Wigner D functions (see eq. (4)in §4. 8.1, [44]) and using the fact that 

D l mm ,{W,p) = e-^d l mml {9)e- m ' p (B5) 
we get the recurrence relations for Wigncr-rf functions: 

sin0d mm , +1 (0j - — - ^ d mm ,(V) (Bo) 

m^(l-m')(l + m ' + l) _ ygj + 1)2 _ gjg _ ml)(l _ m , 

i(i + 1) mm ' 1 J (i + l)(2J + l) mm ' 1 J ' 

Using these relations for d l m2 we may write, 

<4 2 W = MloW + kiO*) + «-i4o(») + «2<Co 2 (f) + «-24o 2 Wl. (B7) 

sin (9) 



where 



m 2 Z 2 -m 2 (Z + l) 2 -m 2 

Kq 



Z 2 (Z + 1) 2 Z 2 (4Z 2 1) (Z + l) 2 (2Z + l)(2Z + 3)' 
y/(l + l) 2 -m? „ V; 2 - m 2 



Ki = 2m — jttttt; Tj k -i = — 2m 



Z(Z + 1)(/ + 2)(2Z + 1)' 1 " " Z(Z 2 — 1)(2Z + 1) ■ 



^[(Z + l) 2 -m 2 ][(Z + 2) 2 -^ 2 ] ; _ v /(Z 2 -m 2 )[(/-l) 2 -m 2 l 



«2 = JZ TTT- — 7-r— — — and K_2 



(Z + 1)(Z + 2)(2Z + 1)(2Z + 3) " ~" 1(1 - 1)(4Z 2 - 1) 

Also, since <4_ 2 (6>) = (- 1 )' +mrf rn 2 ( 7r - and d l m0 (n ~ = {-l) l+m l d l m0 (6) we can write, 

sin (0) 

Using the expression for d m ± 2 we can make the following substitution 

<4 2 (0) + <4- 2 (0) = 2 K MloW + *2^ o 2 (0) + K- 2 4o 2 (0)]/sin 2 0, (B9) 
in the integral we seek to evaluate. We use the following integral for I < I' and L = min{Z, I'} > \m\ > 0, 



d cos 9 



d l m0 (9)d l m0 (9) 



sin 2 ) 



i (i+\ m \y.(i'-\m\y. j 
HV (i-M)!(i'+|m|)! eveni + / 

(BIO) 

k odd I + 1' 



and obviously, for I > I', I and I' have to be interchanged in the above expression. We then obtain 4 + 4-2 
as given in eq. (3.21). 

The integral in eq. (BIO) can also be readily derived. We use the fact that 



d l m o(0) = (-l) m W|^ir(cos0), (Bll) 



which leads to 



f 1 J (9\J' (O\ dcose / (/ — m)\(l' — m)\ f 1 dx 

The symmetry of Associated Legendre Polynomials, -P; m (— x) = (— l) l+m P™ (x) dictates that the integrand 
is antisymmetric for odd values of I + I', hence the integral is zero. However for even values of I + I', we 
can evaluate the integral in the following manner. One of the recurrence relations for Associated Legendre 
Polynomials is ([54], §12.5.) 

Pr{x) = P^ 2 (x) + (21 - l)y/l-x*P™-\x) . (B13) 
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Using equation (B13) we can write, 

i P™ ( x)^(x) T — - 2 = J ^ Pr { x )P^_ 2 { x )-— 2 + (21' - 1) J i P™ ( x)P ; r i 1 ( a; )- 7 =_ . (B14) 

We have provided a proof that the second integral on the right is zero at the end of this section. Thus, from 
eq. (B14) we have 

/! dr I' 1 dr 

i P/"(x)P^(x) T — - = y ^ p r(x) p™ 2 ( X )__^ . (B15) 

In this way we can keep reducing V by two each time until it equals with I (since I + 1' is even and Z < /' it 
will reduce to I). Thus, we have shown that, 



where the second equality follows from the evaluation of a standard integral, which can be obtained, for 
example, from [54]. Substituting in eq. (B12) we can evaluate the integral for m > 0. Clearly eq. (B16) 
is valid for I = I'. For I > I' , I should be replaced by I' in that equation. Moreover, using the property 
d l _ m0 (8) = (—l) m d l m0 (6), we can express the integral for any m ^ 0, as given in eq. (BIO). 

Finally we prove the result used in simplifying eq. (B14) that for even values of I + I' and I < I', 

i Pr{x)P^{x)-f== = 0. (B17) 

Using the recurrence relation of Legendre Polynomials in eq. (B13), we can write 

j^Pr{x)P^ l \x)-^= = J 1 ^P™ 2 (x)P™- 1 Hx)-j^= + (2l-l) J^P^-\x)Pir_- 1 1 (x)dx. (B18) 

Then from the orthogonality relation of associated Legendre Polynomials, 

we can see that the second integral on the right of eq. (B18) vanishes for V ^ I. Thus we have, 

i Pr{x)P^\x)-j== = J J^[x)P^{x)-j=^. (B20) 

We can use the above equation iteratively since the lower indices of the P ; m 's will never match as V > I. So 
the lower index of the first polynomial in the integration can be reduced to either m or m + 1 (depending 
on I — m is even or odd) by repeated use of the above identity. Thus we may write 



,1 , j-l'mW!'-l 

^™ / x ^fin — 1 ' ax « 



j\ P™(x)P™-\x)- 



dx 



Pr(x)P l T_- 1 \x)^== = or (B21) 



y/1 — X 



1-1 C+lW-l 1 ^)^? 



Finally using the relations, 

P% = (-l) m (2m-l)!!(l-x 2 ) m / 2 = (-l)(2m-l)v / l-.T 2 P™r 1 1 , (B22) 
P" + i - x(2m + 1)P™ = x(2m + l)(-l)(2m - 1)^1- x 2 P^ll = (-l)(2m +1)^1- *B23) 

and the orthonormality condition in eq. (B2) we can see that in both the cases right side of eq. (B21) is zero. 
This completes the proof. 
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